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Abstract: A method for designing multiple inputs for real-time dynamic system 
identification in the frequency domain was developed and demonstrated. The designed 
inputs are mutually orthogonal in both the time and frequency domains, with reduced 
peak factors to provide good information content for relatively small amplitude 
excursions. The inputs are designed for selected frequency ranges, and therefore do not 
require a priori models. The experiment design approach was applied to identify linear 
dynamic models for the F- 15 ACTIVE aircraft, which has multiple control effectors. 
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1. INTRODUCTION 

Frequency domain techniques have been used 
successfully to identify dynamic models for aircraft, 
including cases where the aircraft was open-loop 
unstable (Schkolnik, et al ., 1995; Morelli, 2002). 
Recent work has indicated that real-time parameter 
estimation for aircraft dynamic models can be done 
effectively using a recursive chirp-Z transform for a 
selected frequency band, then employing 
equation-error parameter estimation in the frequency 
domain (Morelli, 2000). Advantages of this approach 
include robustness to measurement biases, noise, and 
infrequent dropouts in the time-domain data, 
enhanced signal-to-noise ratio in the frequency 
domain, accurate parameter estimates and error 
bounds in real time, with very low computational and 
memory requirements. 

Good experiment design for dynamic modelling in 
the frequency domain requires that excitation inputs 
to the dynamic system contain a variety of 
frequencies. At the same time, the excitation must be 
such that the amplitudes of the dynamic system 
responses are not too laige, so that the model 
structure assumption, which is typically linear, is not 
violated. For the case of real-time parameter 
estimation on aircraft, inputs are preferably small 
enough so that the pilot cannot distinguish the aircraft 
response to excitation inputs from a typical aircraft 
response to turbulence. Modem aircraft have 
multiple control effectors, so it would be 
advantageous if the excitation for frequency domain 
identification could be applied to multiple control 
effectors simultaneously, so that the amount of time 


required to collect data for dynamic modelling could 
be reduced 

Previous work (Schroeder, 1970) has shown that a 
phase-shifted sum of sinusoids, called the Schroeder 
sweep, provides an input with good frequency content 
and low peak factor. The peak factor is a measure of 
the ratio of maximum input amplitude to input 
energy. Inputs with low peak factors are efficient in 
the sense of providing good frequency content 
without large amplitudes in the time domain. 
Comparisons of the Schroeder sweep with 
conventional linear and logarithmic frequency sweep 
inputs have indicated that the Schroeder sweep is 
generally the superior input for frequency domain 
dynamic model identification (Young and Patton, 
1990). The Schroeder sweep has been used 
successfully in other practical system identification 
problems (Flower, et al ., 1978; Bosworth and Burken, 
1997). 

This paper describes an extension of the Schroeder 
sweep input design method to multiple input design 
with optimised peak factors for real-time parameter 
estimation. The designed inputs are mutually 
orthogonal in both the time domain and the frequency 
domain, and are formulated as perturbation inputs. 
The only a priori information required for the 
multiple input design is an estimate of the 
approximate frequency band for the system dynamics 
and approximate relative control effectiveness for 
proper relative scaling of the input amplitudes. 
Investigations were conducted to determine 
signal-to-noise ratios necessary to achieve good 
parameter estimation results from the real-time 
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parameter estimator. This information is necessary so 
that the input amplitude scaling can adapt to different 
aircraft dynamics throughout the flight envelope. 

The multiple input design technique was applied to a 
lateral / directional linear simulation for the F-15 
ACTIVE aircraft, which includes four control 
effectors. Real-time parameter estimation in the 
frequency domain (Morelli, 2000) was used to 
estimate the dynamic model parameters. Parameter 
estimation results were compiled and analysed. 

The next section describes the multiple input design 
procedure. Following this, real-time parameter 
estimation results from applying the designed inputs 
to the F-15 ACTIVE simulation are presented and 
discussed. 


2 . MULTIPLE INPUT DESIGN 


Each input to the aircraft control surfaces is 
comprised of a set of summed harmonic sinusoids 
with individual phase lags. Each input Uj , applied to 

the 7 th control surface, takes the form 


u j = Z 



( 1 ) 


factors are more desirable for parameter estimation, 
where the objective is to excite the system without 
driving it too far away from the nominal operating 
point, so that model structure assumptions are not 
violated. 

The integers k specifying the frequencies for the 7 th 
input Uj are unique to that input, but the integers k 

used in the summation in Eq. (1) are not necessarily 
consecutive, as will be explained below. The 
objective for the experiment design is to excite the 
aircraft dynamics in a short time period by moving 
multiple control surfaces simultaneously. Since more 
than one surface is being moved, it is advantageous 
for the modelling if the Uj vectors applied to each 

control surface are mutually orthogonal. This helps 
the parameter estimation by completely de-correlating 
the inputs to the aircraft, which improves the 
accuracy of control effectiveness estimates. It is 
possible to make the Uj mutually orthogonal in both 

the time and frequency domains, using inputs 
designed for low relative peak factor, as will be 
shown next. This gives the analyst the flexibility to 
use time domain or frequency domain parameter 
estimation methods, while retaining the desirable 
feature of mutually orthogonal inputs in either 
domain. 


where M is the total number of available 
harmonically-related frequencies, T is the time length 
of the excitation, and the </> k are phase angles to be 

chosen for each of the harmonic components to 
produce a low peak factor PF , defined by 
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where the last equality only holds when Uj oscillates 
symmetrically about zero. In the literature, the 
quantity || Uj || ^ j j| Uj || ^ is called the crest factor. A 

single sinusoidal component from the summation in 
Eq. (1) has = V 2, so the relative peak factor 

RPF, defined by 


RPF(uj) = 


[max^-min^)] 


2>/2 rms^Uj} 


PF( Uj ) 
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(4) 


quantifies the peak factor relative to that of a single 
sinusoid. For a single sinusoid, RPF equals 1. The 
relative peak factor is a measure of efficiency of a 
input for parameter estimation purposes, in terms of 
the amplitude range of the signal divided by a 
measure of the signal energy. Lower relative peak 


In the time domain, a signal composed of a sum of 
sinusoids is orthogonal to any other sum of sinusoids 
with harmonically-related frequencies, regardless of 
the constant phase shift of each sinusoidal component 
contained in the signals. For example, if two inputs 
each contain a single, distinct, harmonically-related 
sinusoidal component, 

(lilt A ( 2 * 4 ‘ A 

«l =«*l —+fa I u 2 =cos T-+& 

then the inner product of these inputs, using the 
discrete-time notation t t -iAt and T ~(N -\)At , is 



= 0 

The sine and cosine of the constant phase angles <j>\ 

and (j) 2 are constants, so the summation equals zero 

because of the orthogonality of harmonically-related 
sinusoids. For more than one sinusoidal component 
in each input, the analysis is similar. So, the inputs 
assembled as in Eq. (1) are orthogonal in the time 
domain. 
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If the frequency indices k that are selected for each 
input Uj are distinct from those chosen for the other 

inputs, then the frequency content of each Uj consists 

of distinct spectral lines in the frequency domain, 
since the summation in Eq. (1) is a cosine series with 
harmonically-related frequencies. Therefore, the 
vectors of Fourier transforms for the inputs as a 
function of frequency have inner products equal to 
zero. In this sense, the inputs are also mutually 
orthogonal in the frequency domain. 

The multiple input design procedure is as follows: 

1. Select the time period T for the excitation, which 
determines the smallest harmonic frequency 
resolution Af = 1 JT and the limit on the minimum 

frequency f mm > 2/7 . 

2. Select the frequency band of the dynamic system 
for the excitation frequencies, [f m j„,f max ] Hz. This 

corresponds to the frequency band where the 
expected dynamic response of the system will occur. 
The frequencies are equally spaced by Af on the 

interval [f min J max ] • The total number of 
frequencies M = fix{(f max - f mi „ )/Af} + 1 , where 

fix indicates rounding to the nearest integer toward 
zero. 

3. Assign approximately an equal number of indices 
k from the set {l,2,...,M} to each input by 

alternating each consecutive frequency among the 
multiple inputs. This approach produces lower 
relative peak factors for the individual inputs, and 
also ensures that each input has frequency content 
distributed evenly across the frequency band 
[ fmm >f max ] • Different assignments of the frequency 

indices could be made for other reasons. For 
example, lower frequency indices might be assigned 
to an input that is known to effectively excite a low 
frequency mode, or particular frequency indices 
might be omitted to avoid exciting an undesirable 
structural response. It was found empirically that if 
the set of selected indices k for a particular input 
consisted of an integer greater than 1, with 2 or 3 
multiples of that integer (e.g., k= 2, 4, 6, or A=5, 10, 
15, 20), the phase angles could be optimised (in 
step 5) so that the relative peak factor for that input 
was very close to 1, and sometimes less than 1. Each 
frequency index can be assigned to only one of the 
inputs, to preserve mutual orthogonality of the inputs 
in both the time and frequency domains. 

4. Generate the input Uj for each of the m controls, 

using Eq. (1) and computing the starting values for 
the phase angles according to method described 

by Schroeder (1970), assuming a uniform power 
spectrum. 


5. Use a simplex optimisation algorithm (Press, et 
ah, 1992) to adjust the <f> k f° r each Uj to achieve 

minimum relative peak factor for that input. The 
optimisation algorithm does not require gradients. 

6. For each input, do a one-dimensional search to 
find a constant time offset for the components of each 
input Uj, so that the input begins and ends at zero 

amplitude. This is equivalent to sliding the input 
along the time axis until a zero crossing is placed at 
the origin of the time axis. The appropriate phase 
shift is added to each sinusoidal component phase 
shift <f> k . Note that to implement a constant time shift 

to all the components, the phase offset for each 
component will be different, because each component 
has a different frequency. Since the components of 
each Uj are harmonics of the base frequency with 

period T, if all the component phase angles <p k are 

shifted so that the initial value of the input is zero, 
then the final value of the input will also be zero. The 
power spectrum, input orthogonality, and relative 
peak factor are all unaffected. 

7. Return to step 5 until the relative peak factor 
reaches a pre-defined goal value or until a maximum 
number of iterations is reached. For this work, the 
relative peak factor goal was set at 1.01 and the 
maximum number of iterations was set to 50. 

Although there are methods for optimising the 
frequency spectrum of inputs for parameter 
estimation, all of them require considerable 
computation, along with some knowledge of the 
system dynamics, usually in terms of a nominal 
model with parameter values. In the current 
application, there is no use for such methods, because 
the intent is to identify the aircraft dynamics in real 
time, for various flight conditions throughout the 
flight envelope, and for arbitraiy failure conditions, to 
enable control reconfiguration. Under these 
circumstances, computing an optimised frequency 
spectrum does not make sense. Instead, the frequency 
spectra for all inputs were defined to be flat across the 
selected frequency band, so that the aircraft dynamics 
would be sufficiently excited, regardless of where the 
modal frequencies happen to be located within the 
frequency band. 

Step 6 of the input design procedure ensures that each 
designed input will be a perturbation, so that the 
designed input can be added to whatever constant 
value the control may have for another reason, e.g., 
trim or manoeuvring. The iteration helps to reduce 
the relative peak factors by adding perturbations to 
the phase angles <f> k (from the phase changes added 
to make the input start and end at zero), which helps 
the optimisation algorithm solve this non-convex 
optimisation problem. 

Figure 1 shows the relative peak factor as a function 
of phase angles for two components in the summation 
of Eq. (1), using 7M5 sec, unit amplitudes, and 
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frequency indices k~2 and &=4. It is clear from the 
figure that this optimisation problem is not convex, 
and also that there are several phase angle solutions 
which are equally good, or nearly so. The data used 
to make the plot in Figure 1 can be used for a global 
exhaustive search for the minimum peak factor in this 
simple case. To generate this data, the relative peak 
factors were computed for values of the phase angles 
on a 2-dimensional grid with intervals of 0.0175 rad 
(1 deg), over the range [ 0, it ] rad ( [ 0, 1 80 ] deg ) . 

The minimum relative peak factor of all these 
computed values was found to be 1.102. Using die 
technique described above to optimise the phase 
angles, the relative peak factor achieved was 1.106. 
This demonstrated that the technique developed here 
found an input with peak factor very close to the 
global optimum, in this simple case. The same 
exercise was then repeated for three harmonic 
components, with frequency indices k= 2, 4, and 6, 
and the same phase angle intervals, in three 
dimensions. In this case, the global exhaustive search 
gave a minimum relative peak factor equal to 1.002, 
while the optimisation technique described above 
produced an input with relative peak factor equal to 
1.003. These investigations suggest that the 
procedure described here does an excellent job of 
designing inputs with very low peak factor. 



Fig. 1. Peak Factors for a Two-Component Input 
3. RESULTS 

For small perturbation motions about a reference 
condition, airplane dynamics can be described by the 
following linear model equations: 

x(t) = Ax(t) + Bu(t) x(0) = * o (5) 

y(t) = Cx(t) + Du(t) (6) 

z i =y i +v i i = 0,l,2,...,N (7) 

Matrices A , 2J, C, and D in Eqs. (5) and (6) contain 
stability and control parameters, which are to be 
estimated from flight data. 

For the linearized lateral/directional dynamics of the 
F-15 ACTIVE aircraft, the state vector x and input 
vector u in Eq. (1) are 


x = [p p r <t>] T u = [S a 8 r 8 dc <J*] r (8) 


The output vector is: 

y = [p P r </> Oyf (9) 

System matrices containing model parameters are: 


( 10 ) 
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Figure 2 shows the NASA F-15 ACTIVE aircraft 
(Doane, et a ! ., 1994). The values of the parameters in 
the model that were used to generate the simulated 
data are given in column 2 of Table 1. 



Fig. 2. F-15 ACTIVE Aircraft 


The multiple input design for the F-15 ACTIVE 
aircraft lateral/directional dynamics includes four 
control effectors: aileron, rudder, differential canard, 
and differential stabilator. Figure 3 shows the time 
histories of the multiple input design for an 18 sec 
manoeuvre, which includes 1 sec for steady trim 
flight at the start, and 2 sec for free response at the 
end. The frequency range for the input design was 
chosen as [0.2, 1.4] Hz. The inputs shown are 
mutually orthogonal in both the time domain and the 
frequency domain. These inputs were designed using 
a flat power spectrum with T=15 sec. Table 2 
contains the amplitudes, frequency indices, and the 
RPF achieved for the 15 sec inputs shown in 
Figure 3. All of the designed inputs achieved an 
extremely low relative peak factor close to 1. 
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Figure 4 shows the simulated responses, which 
include added noise equal to time domain residual 
sequences from flight test data analysis and modelling 
at die same flight condition for the real F-15 ACTIVE 
airplane. This ensures that a realistic noise 
environment was used in the investigations 
concerning the effectiveness of the multiple input 
design. Real-time parameter estimation in the 
frequency domain (Morelli, 2000) was used to 
estimate the parameters. The frequencies chosen for 
the recursive chirp-Z transform were evenly spaced in 
0.01 Hz increments on the interval [0.11, 1.5] Hz. 
This frequency band must include all of the input 
design frequency band, so that the important control 
surface information in the frequency domain is 
included in the parameter estimation. Table 1 
contains the parameter estimation results for this case, 
using only the final values of the real-time parameter 
estimates and standard errors. The parameter 
estimation was started with no a priori information 
for the parameter estimates. The results in Table 1 
show that the model parameters are estimated very 
accurately, and that the computed standard errors are 
representative of the true error in the parameter 
estimates. 

The mean signal-to-noise ratio for the outputs shown 
in Figure 4 was 30. It is of interest to obtain good 
parameter estimation results with the lowest possible 
signal-to-noise ratio on the outputs, so that the aircraft 
and pilot are disturbed as little as possible by the 
activity on the control surfaces. Figure 5 shows the 
effect of output signal-to-noise ratio on the mean of 
the parameter estimate errors. The data for this plot 
was generated by reducing the input amplitudes 
uniformly using the same input forms shown in 
Figure 3, generating new simulated outputs, then 
adding the same noise sequences from flight test data 
used before, so that the output signal-to-noise ratio 
decreased. For output signal-to-noise ratio equal to 
10, the input amplitudes were 0.24 deg, and the mean 
parameter estimate error was 2.7 percent. This choice 
of input amplitudes produced excellent parameter 
estimate accuracy, while keeping the root-mean- 
square value of the lateral acceleration a y below 

0.02 g, which is less than 0.05 g that is typical for 
moderate turbulence. 

4. CONCLUDING REMARKS 

A multiple input design technique for real-time 
frequency domain parameter estimation was 
described and demonstrated. The technique can be 
used to design multiple inputs that are mutually 
orthogonal in both the time and frequency domains, 
with very low peak factors. Input energy over a 
selectable range of frequencies can be efficiently 
injected into the dynamical system by virtue of low 
peak factors and the ability to move inputs 
simultaneously, since the inputs are mutually 
orthogonal. These features make the inputs attractive 
for real-time parameter estimation application, where 
the excitation must have short duration and low 


output response amplitudes, and the dynamic 
response and modal frequencies will change because 
of flight condition changes or failure conditions. 

To use this input design technique effectively in 
flight, the input forms could remain fixed, and a 
simple feedback could be used to scale the input 
amplitudes so that the output signal-to-noise ratios are 
near a chosen acceptable value. The excitation time T 
should be minimized, and this will be impacted by the 
output signal-to-noise ratios. In addition, the F-15 
ACTIVE airplane has a stability augmentation 
system, which will distort the input forms designed 
for parameter estimation. This affects the required 
input amplitudes, and degrades peak factors and input 
orthogonality. Further work is required to address 
these practical considerations. 

The multiple input design technique described here 
can be applied to other problems as well, because 
there is no dependence on a priori models or 
information, other than the expected frequency range 
of the dynamic modes. Specific spectral weighting 
using the component amplitudes can be introduced, if 
warranted by a priori information. 
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Fig. 3. Multiple Input Design 



Fig. 4. Simulated F-15 ACTIVE Response 


Table 1 F-15 ACTIVE Lateral/Directional 



Parameter Estimation Results, 
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Table 2 Multiple Input Design 


Input 

4 

k 

RPF 


0.707 

3, 6, 9, 18 

1.055 

S r 

0.707 

4, 8, 12, 16 

0.995 

Sds 

0.707 

5, 10, 15, 20 

0.995 

Sdc 

0.707 

7, 14,21 

1.003 



Fig. 5. Parameter Estimation Accuracy Dependence 
on Output Signal-to-Noise Ratio 
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